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Abstract — The paper presents an over-parametrization free 
certainty equivalence state feedback backstepping adaptive 
control design method for systems of any relative degree 
with unmatched uncertainties and unknown virtual control 
coefficients. It uses a fast prediction model to estimate the 
unknown parameters, which is independent of the control 
design. It is shown that the system’s input and output tracking 
errors can be systematically decreased by the proper choice 
of the design parameters. The benefits of the approach is 
demonstrated in numerical simulations. 

I. Introduction 

Adaptive control problems are challenging for systems 
with unmatched uncertainties, and backstepping has been 
excessively used to tackle these problems. However over- 
parametrization and ’’explosion of terms” because of repeated 
differentiations are obstacles in its direct applications. One 
way of avoiding over-parameterization is the departure from 
the certainty equivalence principle by modifying the control 
law to include nonlinear damping terms (see [8] for details). 
However, this leads to the adaptation rate to enter into 
the control law and gives rise of high magnitude control 
signals in the case of fast adaptation, which is desirable 
from the perspective of the unknown parameter estimations. 
An alternative way was presented in [2], where certainty 
equivalence control design avoids over-parametrization for 
linear and low relative degree (not exceeding two) nonlinear 
systems in sate feedback settings. 

The ’’explosion of terms” was addressed by combining the 
standard backstepping with a sliding mode control approach 
[1], [4], by the Multiple surface sliding control [17], or by 
the Dynamic scarface control [16]. Alternative approaches 
include approximation of the virtual control derivatives using 
sliding mode filters [9], [15], neural networks [11], [12], 
fuzzy systems [5], first order linear filters [18], and second 
order linear filters [3], The authors of [3] used the singular 
perturbation method and Tikhonov’s theorem to prove the 
closed-loop stability and to obtain the performance bounds. 

In [14], we introduced a certainty equivalence state 
feedback indirect adaptive control approach without over 
parametrization for nonlinear systems of any relative degree 
in the parametric strict feedback form. The approach was 
based on the state prediction model, which is capable of 
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providing fast estimation of unknown parameters indepen- 
dent of the control design. This property is the consequence 
of feeding back an error term with the gain proportional to 
the square root of the adaptation rate, like in the modified 
reference model MRAC (M-MRAC) architecture introduced 
in [13], where it has been shown that the error feedback gain 
acts as a damping factor for the adaptive signals, while the 
adaptation rate determines their frequency. 

In this paper we extend the method to the systems with 
unknown virtual control coefficients. It is shown that the 
input tracking error (difference between ideal control and 
command filtered certainty equivalence control signal) and 
output tracking error can be regulated as desired by the 
proper choice of design parameters. 

The rest of the paper is organized as follows. In Section II, 
we give the problem statement and the main assumptions. In 
Section III, we introduce the identification model and give its 
properties. The control design and performance analysis are 
presented in Section IV. A simulation example is presented 
in Section VII, and some concluding remarks are given in 
Section VIII. 

II. Problem Statement 

We consider control design for a system 

±i(t) = ctiXi + i(t) + & J 

i = 1, . . . ,n - 1 (1) 

X„(t) = a n u(t) + 'dlf n {t ) , 

with some initial conditions tc(0) = Xo, where x £ R n and 
u £ K. are the state and input of the system, $ £ is 
a vector of unknown constant parameters, f l : K. — > 
i = 1 are sufficiently smooth known functions, and 

a±, ... ,a n are unknown nonzero virtual control coefficients. 
We assume that the sign of a n is known (positive without 
loss of generality). Here, we use short hand notations f ^t) = 
■ ,Xi(t)), i = l,...,n. 

Our goal is to design a controller such that all closed- 
loop signals are bounded and the output y(t) = X\ (t) of 
the system (1) tracks a desired bounded command, which 
has bounded derivatives up to of order n — 1. One way 
of providing such a command is to filter appropriately cho- 
sen piecewise continuous signal y com (t) trough a reference 
model of the same order as the system itself, and track the 
model’s output, which is the approach adopted here. The 
reference model is given in the controllable canonical form 

= kl rn,x m {t) -T b m y corn (t) , x rn {0') = Xq , (2) 



The prediction error dynamics are easily derived to be 


where x m £ R" is the state of the reference model, A rn = 
A v n ^ n k x , A — [0 v n i . . . v n , n —i], b m — k r v nn , 
v Ut i is the i-th coordinate vector in R", and the gains k x 
and k r are chosen from the perspective of the performance 
specifications and to make A m Hurwitz. The output of the 
reference model is y m (t) = x m p(t). 

III. Identification 
We transform (1) into the form 

±i{t) = x i+1 (t) +biX i+1 (t) 

i = 1, . . . , n — 1 (3) 

x n (t) = g n (t ) + a n [u(t) - b n g n (t) + dff n (t )] , 

where bi = a*- 1, i = 1 , . . . , n- 1, b„ = l/a„, d„ = b n d n , 
and g n (t) is a smooth function to be specified in the control 
design. This form is better suited for the parameter estimation 
scheme adopted here and does not require devision by the 
parameter estimates in the control design. 

To estimate the unknown quantities in the first n — 1 
equation in (3) we design an identification model 

xft) = x i+1 (t) + bi(t)x i+1 (t) + + Xxi(t) 

i = l,...,n-l, (4) 

where xft) is the prediction of the z-th state, xft ) = 
xft) — xft) is the prediction error, A > 0 is a design 
parameter, bft), dft) are the estimates of the unknown 
constant parameters bi, r>, . These estimates are generated 
according to adaptive laws 

bi(t) = 'yx i (t)x i+1 (t) (5) 

&i{t) = 'yx i (i)f i {t), i = l,...,n-l, 

where 7 > 0 is the adaptation rate. 

The identification model for the estimation of the uncer- 
tainties in the last equation in (3) has the form 

Xn{t) = g n (t) + \x n (t) (6) 

+ a n {t) u{t) -b n (t)g n (t) +&1{t)f n (t) , 

where x n (t) is the prediction of x n (t), x n (t) = x n (t)—x n (t) 
is the prediction error, a n (t), b n (t), are the estimates 

of the unknown constant parameters a n , b n , d n . We notice 
that if we design the control signal as 

u(t) = b n (t)g n (t) - tfl(t)f n (t) , (7) 

then the prediction model (6) reduces to 

Xn(t) = g n (t) + A x n (t) . (8) 

Therefore a n (t) is not needed for the control design, and for 
the generation of the necessary estimates b n (t) and i9 n (f), 
the adaptive laws 

b n (t) = -7 x n (t)g n {t) 

&n(t) = 7 x n (t)f n {t). (9) 

can be implemented with x n (t) generated using the reduced 
model (8). 


x(t) = (A — Xl)x(t) + r)(t) , (10) 

where %(f) =b i (t)x i+ i(t) + : d i (t)f i (t), i = l,...,n-l, 
Vn{t) = a n [—b n (t)g n (t) +d n ( t )fn ( *)], and the parameter 
estimation errors are defined as bi(t) = bi — bi(t), i = 

1 . . . ,n, = ■di - &i(t), i = 1 ... ,n - 1, tf n (t) = 

tin-Vnit). 

Lemma 3.1: The error signals Xi(t), i = 

1 ,n are globally bounded, and Xi(t) — > 0 as t — > 00. 
Proof: Consider a candidate Lyapunov function 

n 1 n— 1 

y {t) = '52 i i( t ) + -'52 

i= 1 ' i— 1 L 

7 L J 

It is straightforward to compute the derivative of V (t) along 
the trajectories of the prediction error dynamics (10), and the 
adaptive laws (5) and (9) 

n n 

V(t) = -2A ^2xi(t) + 2y^x i - 1 (t)x i (t) . (11) 

i = 1 i=2 

Since 2x i _ 1 (t)x i (t) = [zj_i(f) +Xi(t)] 2 - [xf^t) +x? (f)], 

we conclude that 

n 

V(t) = -{2\-l)^x?{t)-[xl(t)+xl(t)\ (12) 

i = 1 

n— 1 n 

-J2[xi(t) -x i+l (t)] 2 < -(2A- 1)^2 x?{t). 

i= 1 i—1 

If we select 2A > 1, it follows from the LaSalle-Yoshizawa 
theorem ([8], p.24) that x(t). Gift), bi(t), i = 1, . . . ,n are 
globally uniformly bounded, and x(t) — > 0 as t — > 00. ■ 

When x(t) and u(t) are bounded (which will be provided 
by the control design), the following upper bounds on the 
error signals x(t) and rift), i = 1 , . . . , n can be derived. 
Lemma 3.2: Let the estimates xft), bft), and ‘Oft), i = 

1. . . . , n be generated by the identification models (4), (5), 
(8), and (9). In addition, let x(t) and u(t) be bounded. Then 
rjft) and xft) satisfy the following bounds 

\Vi{t)\ < Pi,\e- Vt + (13) 

\xi{t)\ < Pi^e- Vt + ^ , (14) 

where the constants ffj > 0, i = 1 j = 1,2, 3,4 

and v > 0 are defined in the proof. 

Proof: Differentiating fg(t), i = 1, ...,n — 1 and 

substituting the adaptive laws we obtain 

Vi(t) = bi(t)x i+1 {t) +-dj(t)fi(t) + bi(t)x i+1 (t) 

+ ) = -7 Pi(t)xi(t) + hft ) , (15) 

where pft ) = xf +1 (t) + fj ( t)f ft ), hft) = bft)x l+1 {t) + 
■di ( f)fi(t ), i = 1, . . . , n — 1. It is straightforward to show 



that the second derivative of the signals 77* (t) , i = 1 , . . . , n— 
1 satisfy the dynamics 

m(t) + + 7 Pi(t)Vi{t) = (16) 

- 7 Pi(t)Xi(t) + A hi(t) + h i(t) . 

Similarly, for rj n (t) we obtain 

rin(t) + Xfln(t) + 7 Pn(t)rin(t) = (17) 

- 7 Pn(t)x n (t) + A h n (t) + h n (t) , 

where p n (t) = a^gl{t) + fZ(t)f n (t)\ and h„(t ) = 
a n [-b n (t)g n {t) + (t)f n (t)]. 

Since x(t) and u(t) are bounded, and g(t), i = 

1 , ,n are smooth functions, there exist positive constants 
Si, 1, Si , 2 , Si , 3 such that ||p i (f)||,c 00 < 5 it i, ||p i (t)||,c 00 < 
Si , 2 and ||/ii(f)|| / c oo < c/3. Then, it follows from the 
results of [13] that choosing A > 2^/Soj, where So = 
max^i, . . . , S Ui 1], damps the oscillations in f]i(t) and guar- 
antees the bound 


\Vi(t ) I < fte 


— t 


AaPWII + ^IM*)I , 


y(f) < — (2A - 1 ) 

Integrating (19) we conclude that 

V{t)< 


V(t) - -Si 
7 . 


F(0) - f 


e -(2A-l)t + 

7 ' 


Since ||ce(t)|| 2 = J2i= 1 *?(£) < V(i), we obtain 

ll*(*)ll < 


S 2 

21/(0) - -7- 

7. 


5 2 


e -(2A-l)t + 71 > 

7 


|a(i)|| < 5 2 e- (A ^ 1/2)t + ^ 


a/7 


®n(<) = -\x n (t) + T] n (t) . 

Pn, 


Therefore 

\x n {t)\ < \x n (t)\e~ xt 

Pn, 2 


A- 


1 r 

I 


e _!/t - e -At l 


A^7 


[1 - e~ xt ] < p n , 3 e- vt + 


where and /3 Hi 4 are readily computed. Working back- 
ward and taking into account the bounds 


\x i+1 (t) + rji(t)\ < (Pi , 1 + Pi + i,3) e 

we obtain in the similar way the inequalities 

„ — vt 1 PiA 


-ft A Pi t 2 + pi+ 1,4 


A^7 


\xi(t)\ < Pi, 3 e~ 


Av^ 


(25) 


(18) 


where v\ is proportional to y/y, and the positive constants 
Pi and S iA are independent of 7 (see details in [13]). 

It follows from Lemma 3.1 that there exists 5\ > 0 such 
that ]C”=i[ll^iWI| 2 + Pf(t)\ < Sf. Therefore, the inequality 
(12) can be written as 


(19) 


( 20 ) 


( 21 ) 


Recalling that y/a + b < y/a + a fb for any a > 0, b > 0, we 
arrive at 


( 22 ) 


where S 2 = yj 1 21/(0) - ^|. 

Substituting (22) into (18) we obtain (13) with v = 

min(i/i, A - 1/2), j3 iA = Pi + S it2 S 2 and = Si, 2 S\ + 
Si, 3 Si,4’ 

The bound (14) is obtained by direct integration of (10), 
the last equation of which has the form 


(23) 


(24) 


for each i = n — 1, . . . , 1. The proof is complete. ■ 

IV. Control Design for Known Systems 

In this section we design two controllers, which are solely 
used for the performance analysis. 

A. Conventional Backstepping 

We first design a controller assuming that b, 
are known (ideal control) and formally applying the conven- 
tional backstepping procedure to the system (3). The error 
variables are defined as 

e?(f) = x°(t) - x mA (t) (26) 

e i(t) =x°i(t) -x m ,i(t) -g^it), i = 2 , . . . ,n , 

where the stabilizing functions have the form 

9 i(t) = -bix° 2 (t) - 

9i(t) = — bix° i+ 1 (t ) - tij f°(t) + g^^t), (27) 

i = 2 , . . . , n — 1 . 

To distinguish between the designs we furnish all variables 
with a superscript 0, which indicates that the loop is closed 
by the backstepping controlled designed for (3) in the ideal 
case of known dynamics. 

The eiror variables satisfy the equations 

e°i(t) = e° +1 (f), i = 1 , . . . ,n — 1 (28) 

e° n (t) = gl(t) + [u\t)-b n gl(t) + tiZf n (t)] 

- Xm,n(t) ~ g n -l(t) ■ 

When the control signal is designed as 

u°(t) = b n g°(t) - &Z f n (t) , (29) 

the error dynamics reduce to an exponentially stable system 

e°(t) = A m e°(t) , (30) 

if we select 


9n(t) = ~K e (*) + Xm,n(t) + 9n-\(t) ■ 


(31) 


A\/7 ’ 


Although the stabilizing functions depend on the corre- 
sponding virtual controls, it can be easily shown that the 
algebraic loops are solvable and the controller could be 
implemented if the uncertainties were known. 

We notice that since e°(t) is exponentially stable, it 
follows that y°(t) exponentially converges to y m (t). In 
addition, y CO m(t) €E £00 implies that x m (t) £ £00, which 
along with e°(f) € £00 can be used to recursively show that 
5°(f), 9i(t) G £00, i = 2, ...,n. Then, u°(t) G C x , and 
the following lemma has been proven. 



Lemma 4.1: The controller defined by (26), (27) and (29) 
guarantees the boundedness of all closed loop signals in the 
system (3) and exponential tracking of the reference model’s 
output. 


B. Command Filtering 

Following [3], we introduce the command filtered version 
of the design (26), (27) and (29). In this case the error 
variables are introduced as 

e{ (t) = x{ (■ t ) - x m ,i (t) (32) 

4 (t) = x{ (<) - - of_ 1;1 (f), i = 2,...,n, 

where (t), i = 1, ...,n — 1 are the outputs of the 
command filters 

4,1 W=wcr° 2 (f) (33) 

4,2(1) = -2Cwcr ° 2 (f) - ^>[4,i(t) - g{ (f)] 

with initial conditions 0 ® x (0) = gf (0) and of 2 (0) = 0. The 
stabilizing functions are defined similar to the ideal case as 

gl(t) = -bix£(t)-‘&Jf[(t) (34) 

4 (' * ) = -hx{ +1 (t) - f{ ( t ) + WCT°_ lj2 (f), 
i = 2, . . . , n — 1 . 


where the superscript / indicates that the corresponding 
quantities are computed with the command filtered control 
in the loop. The system (3) is written in error variable as 


4(t) = 4+i(t)-4(t) + 4i(t), i = 1, ... ,rz. — 1 

4(t) = 9n(t) + a n u(t) - b n g f n (t) + 

Xm,n(t) ^^n—lgit) • (35) 

The command filtered control signal is defined according 
to equations 


u f (t) = b n g S n (t) - •&1 f{ (t) (36) 

ffn(t) = -kje f (t) + x m , n (t) + u;o 0 n - 12 (t) , 

which translates the error dynamics (35) into 

e f (t) = A m e f {t) + a(t ) , (37) 


where 


a(t) 


4,i(t)-9i(t) 

4-i,i(t) - 4-i(t) 

0 


For the stability analysis we introduce the compensated error 
signal as e{(f) = e-'(f) — £(t), where £(t) satisfies the 
dynamics 


Obviously, the compensated error dynamics are exponen- 
tially stable. 

Lemma 4.2: The command filtered controller defined by 
(32), (33), (34), and (36) guarantees the following relation- 
ships 

e f(t)-e 0 (t) = O(e), m=0(e) (40) 

4,i(t) - 9i(t) = O (e) , i = l,...,n-l 
u4,2(t) ~ 9i(t) = ° ( e ) > 

where e = 1 /oj (the proper choice of ( and u> is discussed 
in [3]), and the notation O (e) is adopted from [6] (p. 383). 

The proof of the lemma follows the steps of Theorem 2 
from [3]. 

V. Adaptive Control Design 
A. Certainty Equivalent Control 

Here, we design a controller for the identification models 
(4) and (6). In essence, the form of the control signal is given 
by (7) and we need to only specify g(t). To distinguish from 
the ideal design we use ”hat” notation for the corresponding 
variables. Following the conventional backstepping design 
steps we define the stabilizing functions as 

9i(t) = ~bi(t)x 2 (t) - (t)f 1 (t) 

9i(t) = -k{t)x i+ i(t) - (t)fi(t) + gi-i(t ) , 
i = 2, , , . , n — 1 . (41) 

and the error variables as 


ei(t) = xi(t) - x m> i(t) (42) 

e,(f) = Xi(t) - x m ,i(t) - gi-i(t), i = 2,...,n. 

The identification model (4) in the error variables takes the 
form 


i(t ) (f ) T (f), i 2 , ... ,n 1. 

The dynamics of e n (t) is derived from the equation (6) as 
follows 

e„(t) = g n (t) + a n u(t) - b n {t)g n (t) + &n(t)f n (t) 

- X m ,n(t) ~ 9n-l(t) + \x n (t) . (43) 

We notice that the error system is translated into 


e n (t) = —kje(t) + Xxi(t) 


(44) 


which is asymptotically stable, if we set g n (t) = — kj e(t) + 
Xm,n(t) + g n -i(t) and the certainty equivalence control law 

u(t) = Kit)g n {t) - 'd T n ( t )f n ( t ) ■ ( 45 ) 

It is easy to see that the error dynamics can be written as 
e{t) = A m (t)e(t) + Xx(t) . (46) 


&(t) = Anti(t) + a(t) (38) 

with the initial condition £(0) = 0 and £ n (t) = 0 for all 
t > 0. It is straightforward to verify that 


Lemma 5.1: The controller defined by (41), (42) and (45) 
guarantees boundedness of all closed-loop signals and the 
inequality 


Pi, i 


v 2 X^ ' 


4(t) =A m e{(t). 


(39) 


\e t .(t)\ < p ij5 e V2t 


(47) 



where /3jj > 0, i = 1, . . . , n, j = 5, 6 and j/ 2 > 0 are 

defined in the proof. 

Proof: According to Lemma 3.2, x(t) £ Coo. Since A rn 
is Hurwitz, it follows from (46) that eft) £ Coo ■ In addition, 
from y CO m(t) £ Coo it follows that x m (f) £ Coo- Therefore 
xi (t) = ei (f) + x m i(t) £ Coo and x\ (t) = Xi(t) + x\ (f) £ 
Coo- Then gi(t) £ Coo, since fi(xi(t)) and hi(xi(t)) are 
continuous and 6i(f), 9i(t) £ Coo according to Lemma 3.1. 
This implies that x 2 (t) = e 2 (t) + g\(t) G Coo and ;r 2 (f) = 
a; 2 (t) + i 2 (f) £ Coo- Continuing this recursion we conclude 
that x(t) £ Coo, x{t) £ Coo, 9 i{t) £ Coo, i = 1 ,---,n, 
gft) £ Coo, i = 1, ■ • • ,n — 1, and u(t) £ Coo- Thus all 
closed-loop signals are bounded. To derive the bound (47) we 
notice that e Amt < ne~ U2i for some k > 0, where ^ 2 > 0 is 
the decay rate of the reference model. Therefore integrating 
(46) and taking into account (14) we obtain 


|e*0)| < 


K0j, 3 
V — V 2 


e — e 


,-^ 2 1 _ -vt ] 






-V-2 1 


Pi, 6 
vixs/a 


[l - e~ V2t \ 


(48) 


This control design is used for the performance analysis. ■ 


B. Actual Control Design 

Now we design the actual control applied to the system 
as the command filtered version of the certainty equivalent 
control from the previous subsection. The uncompensated 
error variables are introduced as 

e{ ( t ) = x{ ( t ) - x mi i (f) (49) 

e{(t) = Xi(t) - x m ,i(t) - <7j — i = 2, . . . ,n. 

where is the first state of the command filter 


reduces the error dynamic to 

e f {t) = A m e f (f) + a(t ) + A x f (t) , (54) 

where a(t) is obtained from ot(t) by replacing of jff) with 
07 , 1 0 ) and g{(t) with g{(t). 

The compensated error signal is e{,(t) = e^(t) — £(f), 
where £(f) now satisfies the dynamics 

(ft) = A m ift) + a(t) (55) 

with the initial condition £(0) = 0. It is straightforward to 
verify that 

e{ (f ) = A m e f c (f ) + Xx f ( t ) . (56) 

Lemma 5.2: The controller defined by (49), (50), (51), and 
(53) guarantees the following relationships 

e f (t)-e(t)=0(e), £(t)=0(e) (57) 

(t) -g{(t)=0(s), i = l,...,n-l 

w(j j, 2 (f) - g{(t) = O (e) , i = 1, . . . ,n - 1 . 

Proof: Since the exponential convergence of xJ (t) 

is not guaranteed, Tikhonov’s theorem ([6], Theorem 11.2) 
cannot be directly applied to the system comprised of (52), 
(55) and (56). However, since x* (f) does not depend on e, 
the state transformation s^(t) = (f) — n(t), where /i(f) 

is dynamically defined as /7(f) = A m /r(f) + Xx^ (f), results 
in the system 

s f (t) = A m s f (t) + a(t) . (58) 

Next, define s(t) = e(f) — fi(t), which satisfies the expo- 
nentially stable dynamics 


07 , 1 (f) = w<r i>2 (t) (50) 

2 (f) = - 2 C £*> 07,2 (i) - w[<J,:,i(f) - g{ (f)] 

i = 0, . . . , n — 1 , 

with the initial conditions (7,,i(0) = g{ (0) and <7,, 2 (0) = 0. 
The stabilizing functions have the form 

9i(t) = -bi{t)x{{f) -df(t)f{{t) 

9i (*) = ~k(t)x{ +1 (t) - ( t)f{(t ) + WCr,_i j2 (f) , 

i = 2, . . . ,n — 1 . (51) 

where the superscript / indicates that the corresponding 
quantities are computed with the command filtered control 
in the loop. The identification model in error variables takes 
the form 

ef (f) = e i+1 (t) - g{ (f) + ^^(f) + Xx{ (f), 

i = 1, . . . , n — 1 (52) 

el(t) = gl(t) + [ u(t ) - b n (t)gl(t ) + 

(f) + UJO n - 1, 2 (f) + Xx^(t) . 

where x^ (t) = X* (t) — X* (f). Designing the controller as 

u f (t) = b n (t)g{{t)-V T n (t)f f n (t) (53) 

gl(t) = -fcfe / (f) +x m ,n(t) + ■ 


s(t) = A m s(t) . (59) 

Following the steps from [3], it can be verified that 
Thikhonov’s theorem’s conditions are satisfied for the sys- 
tems (58) and (59). Therefore the last three relationships in 
(57) can be concluded along with s^(i) — s(t) = 0(e). 
Then, it follows that e^(f) — e(f) = s^(t) — s(t) = O (e), 

which completes the proof. ■ 

VI. Performance Analysis 

Lemma 6.1: Let the command filtered controller for sys- 
tem (3) be defined by (32), (33), (34), and (36). Then all 
closed-loop signals are bounded and 

x^(t) — x°(t) = O (e) (60) 

u^(t) — u°(t) = O (s) . 

Proof: From Lemma 4.1 we have e°(f) £ Coo, 

9i(t), 9i(t) G Coo for i = 1. Therefore, (40) 

implies that e* (t) G Coo, of (t) G Coo an d <rf 2 (i) G Coo 
for i = 1,. . . ,n — 1. Next, (32) implies that x * (t) £ Coo- 
Therefore f{(t) £ Coo, an d equations (34) imply that 
g{ (f) G Coo for all i = l,...,n. It follows from the 
definition (36) that w(t) £ Coo- Further, we notice that 

x{ (' t ) - x°(t) = e{ (f) - e°(f) + erf i(f) - g- (f) . 



Since e{(t) - e°(t) = O (e) and a -^(i) - g?(t) = O (e), it 
follows that x{(t) — Xi(t) = O (s) for all* = 1, . . . , n. From 
Lemma 4.2 we also have ujo n -i, 2 (f) — 9n-i(t) = ® ( £ )> 

therefore 

9n(t)-9n(t) =°( £ ) ■ 

It follows from the smoothness of f n that f n (x* (f)) — 

fn( x °(t)) =c, ( £ )- Then 

— u°(t ) = O (e) . 

The proof is complete. ■ 

Lemma 6.2: Let the controller for system (3) and the 
identification models (4) and (8) be defined by (49), (50), 
(51), and (53). Then all signals are bounded and 

x^(t) — x(t) = O (e) (61) 

— u(t) = O (e) . 

Proof: From Lemma 5.1 we have e(t) £ Coo, 

gift), gi(t) £ Coo f°r i = 1 , ,n — 1. Therefore, Lemma 
5.2 implies that e-f(f) € Coo , G Coo and (Ji.-iit) £ 

Coo for i = 1, . . . , n — 1. It follows from (49) that x^(t) £ 
Coo ■ Since (t) £ Coo, it follows that x? (t) £ Coo as 
well. Therefore (t) £ Coo, and equations (51) imply that 
g{ (t) £ Coo for all* = 1, , n. The definition (53) implies 
that u-f(f) £ Coo as well. Next, we observe that 

x{ (f) - Xi(t) = e{ (f) - ii(t) + - g.i{t) . 

Since e{(t) — Sift) = O (e) and (Ti,i(t) — gpt) = O (e), 
it follows that x { ( t ) — xfi) = O (e) for all i = 1, ... ,n. 
On the other hand, wovi- 1 . 2 (f) — 9 n-i{t) = O (e) is true 
according to Lemma 5.2, therefore 

i/nW ~9n{t) =0{e) . 


Next, f n (xf(t)) — f n (x(t )) = O (e) holds because f n is 
smooth. Since b n (t) and i? n (f) are bounded, it follows that 
u-t(f) — u[t) = O (e). The proof is complete. ■ 

Now we are ready to prove the main result. 

Theorem 6.1: Let the controller for the system (1) be 
defined according to command filtered scheme given by (4), 
(5), (8), (9), (49), (50), (51), and (53). Then the input and 
output tracking errors satisfy the following upper bounds 


I y(t) - Vm(t)\ < fae V2t - 
|fi'(t)-u°(t)| <fae~ V2t 


fa 


+ 0(e) 



+ O (e) , 


(62) 

(63) 


where @ 3 , pi, @ 5 , /3q, and u 2 are positive constants defined 
in the proof. 

Proof: First of all, we notice that y(t) — y m (t) = 

x{(t)-x ml (t) = x({t) + x[{t) - x ml (t) = x{(t) + e{(t). 
Since S{(t) = e± (t) + O (e), it follows that y(t) — y m (f) = 
x{(t) + Si (t) + O (e). Using (14) and (47) one can obtain 



since v > v 2 for fast adaptation, and @3 and Pi are readily 
computed. 

To compute the bound for the control signal, we consider 
two closed-loop systems. Namely, the system (1) with the 
controller defined by (32), (33), (34) and (36), and the system 
(1) with the actual controller. For the clarity, we denote 
y(uf), when y(t) is generated by the controller u^(t), and 
y(u^) when y(t) is generated by the actual control u^(t). 
The error between this variables is yf(t) = y(u^) — y(u f). 
First, we recursively compute x{(t) = xfu*) — xfiJ), 

fi(t ) = fi(u f ) ~ fi(u f ), g{ (t) = 9 i{u f ) - 9 i(u f ), i = 

1 For i = 1 we have x{(t) = Xi(iC) — x m> i (t) + 

x m} 1 (t)-x 1 (u f )+x 1 (u f )-xi(u f ) = e{(t)—e{(t)—xi(t) = 
e°(f) + 0 (e) — xi(t) — e{(f). Since |e°(f)| < |e°(0)|Ke _iy2 *, 

it follows from (14) and (47) that 

\x{ (f)| < || e° ( 0 ) || Ace - + 0 (e)+ A, 3 e _iy * + 
+fa,i>e- v * t + ^ < Px,i e ~ V2t + ^ + O (e) , (64) 

where f3 x>1 = || e° ( 0 ) || « + /?i , 3 + /3i >5 and S x>2 = + 

Since x(f ) is bounded and is smooth, it follows from 

(64) that 

ll/i (f)ll — fa,i e ~ U2t + ~£= + O (e) . (65) 

V7 

where (if a = LiP x \, Sf t i = LiS X 2 , and Li is the Lipschitz 
constant for f l . Using (34) and (51) it is straightforward to 
compute that 

9i(t) = ~bix{ (t) - - 91 (f) . ( 66 ) 

On the other hand, it follows from (32) and (49) that 

x -2 (t) = e{(t) + <Ti,i (u s ) - x 2 (t) - e{(t) - 01 ,i(u f ) ■ (67) 

From Lemmas 4.2 and 5.2 it follows that eri i(iU) = 
gi(uf) + O (e) and <Ji,i(uf) = gi(u^) + O (e), therefore 

x{ (f) = 4 ( t ) - x 2 (t) - e{ (■ t ) + O (s) 

-bix{(t) - ‘&lf[(t) - r)i(t) . ( 68 ) 

Solving ( 68 ) for xf 2 (t) and recalling that 1 + 61 = aiO, we 
conclude 

4 (*) - £ 2 ( t ) - e{ (t ) + O (e) - &J f{ (t) - 771 (t ) 
x 2 \ l ) — ’ 

a 1 

from which we obtain 

l^(*)l < \k\ [l|e°( 0 )||«e— 2 * + fa, 3 e-* + 

+fa,Be~’* t + 4 ^ + ll^i II [/3/,ie-" 2t + (69) 

+fa,ie-* + 2 %) + O (e) < p x , 2 e-^ + ^ + O (e) , 

where fi x / 2 and 5 Xj2 are readily computed. Taking into 
account (69), we conclude from ( 66 ) that 

\gM < \bi\ [fate-** + + lltfill [p f ,ie~ V2t + 4 + 

+/3i,ie— * + ^ + O (e) < P g ,ie~^ + + O (e) . 



To prepare the next step, we observe that the signals 

crj,i(t) = of iO) - 071(f) and CT», 2 (t) = ct° 2 (*) - 072(f) 

satisfy the operator equations 

^,i(s) = Gi(s)s/ (s), Gi(s) = a2+2 ^ fl+tja (70) 
ct 4 , 2 (s) = G 2 {s)g{ (s), G 2 (s) = ,, 2 +2 ^ sW • 

Since ||Gi (s) = 1 for ( > V2/2 and ||G 2 (s)||« 00 = 
(2£) _1 , it follows from (70) that 

|<7i, 2 (f)| < Ar,l e 1/24 4 ~h + O (e) , (71) 

where /3 CTj i = (2C) -1 /? s ,i and <5 CTi1 = (2C) _1 5 ga . Using 
again (34) and (51) we compute 

92 (*) = -hx f 3 (t) - i?J/ 2 (f) - 772(f) + wcr 1)2 (f) . 

Now, we use (32) and (49) to write 

x{ (f) = e{{t) + cr 2) i(u / ) - x 3 (t) - e 3 (t) - ct 2)1 (u / ) , (72) 

where (J 21 (u^) — cr 2il (u^) = g{(t) + O (e) according to 
Lemmas 4.2 and 5.2. Therefore 

4 ( t ) = e{ (f) - x 3 (f) - 4 ( t ) - b 2 x f 3 (f) - f f 2 ( t ) 

— 772 (f) + W(Ji, 2 (i) + O (e) , (73) 

which after solving for x^f) reduces to 

(f ) = ^ [4 (f ) - ®3 (f ) - 4 (f ) - i? 2 fl (; f ) - t? 2 (f ) 

+w<7i >2 (f) + 0(e)] , (74) 

~ f 

wher e/ 2 (f) satisfies the bound 

||/ 2 (i)|| < Pf,2 e l ' 2t + -^ + O (e) , (75) 

VT 

which follows from the smoothness of / 2 (-) with respect to 
x \ and x 2 and from the bounds (64) and (69). Substituting 
the corresponding bounds in (74) we obtain 

\4(t)\ — + —j= + 0(e) , (76) 

which enables us to compute the following bound for gi, (t) 

14(f) I < N [p x , 3 e-^ + +/3 2 ,ie-^ + ^ 

+ ||« 2 || + ^] + t [P^ e ~ U2t + 


°x,3 , /O l/t P2,2 

77V +P2.ie +-7= 


p-^2* I £/ 


\/l \ + 2C 


^ ' U7 

— i/?t I G, 1 

V7. 


+0 (e) < P g , 2 e- V2t + 5 -ff+0(e). (77) 

Continuing in the same manner, we can derive the follow- 
ing bounds 

II /f toll < ^‘ + t + 0 ( £ ) 

|®f(t)l < Px,ie-~ U2t + + o (e) 


l4tol < /3 s .ie 1224 + -^= + G (e) 


for all* = 1 , ... , n. 


Now, we are ready to derive a bound for the control signal. 
To this end, we recall that u^(t) — u°(t) = O (e) according 
to Lemma 6.1. Therefore 

u°(t) — u^(t) = 0(e ) + v,f (f) — v,f(t) . (79) 

Substituting the expressions of w(t) from (36) and u'(i) 
from (53) into (79) we obtain 

u°(t) - u f (t) = O(s) - b n g n {u f ) - &Zf n {u f ) 

+b n (t)g n (u f ) + fi T n (t)f n {u f ) (80) 

= ~Vn(t) ~ b n g n (t) - &Zfi(t) + 0{e), 

which upon substitution of the corresponding bounds results 
in 

I U°(t) - u f (t)\ < \b n \ /3g, n e- V2t + + h,ie~ vt 

+^ + IKH [/3/,„e-^‘ + *-**] +0(e) 

= + ^ + 0(e), 

which completes the proof. ■ 

Remark 6.1: It follows form Theorem 6.1 that the bounds 
on the input and output tracking errors can be systematically 
decreased by choosing large values for u> and 7. It should 
be notice that u> appears in the bounds of the control signal 
not only in O (e) term as e = — , but also in the coefficient 
/?6 in the form of a factor w n . Therefore in order to make 
the term smaller the adaptation must be much faster than 
the command filter’s response. 

VII. Numerical Simulations 
As a simulation example we consider a third order system 
with fx(xi) = 2.2x1, f 2 i. x i > x 2 ) = 1.3xix 2 and f 3 (x) = 
0. For the numerical simulations the uncertain parameters 
are set to -t?i = 2.2, r d 2 = 1.3, a± = 1.2, a 2 = 2.4 and 
a 3 = 2.5. The reference model is implemented with k x = 
[8 10.4 5.2] t and k r = 16. The frequency and damping ratio 
of the command filters are respectively set to w = 500 and 
C = 0.8. The adaptation rate is set to 7 = 50000 and the 
error feedback gain is A = 2^/7. 

Figure 1 displays the controller performance for the exter- 
nal step command of magnitude 2. Clearly asymptotic track- 
ing is achieved for both output and input signals with a small 
transient error. System’s sinusoidal response is presented in 
Figure 2. Good tracking performance can be observed in this 
case as well. 

VIII. Concluding remarks 
We have presented a certainty equivalence adaptive back- 
stepping control method for nonlinear systems with un- 
matched uncertainties and unknown virtual control coeffi- 
cients without over-parametrization. The approach uses a fast 
identification model, which is independent of the control 
design, and the command filtered backstepping method. 
This separation of the estimation and control design is an 
important feature that enables the designer to achieve desired 
transient and steady state properties by the proper choice of 
the control parameters. 
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Fig. 1. Step response and the corresponding control signal history. 
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